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Abstract: We study the M = 4 theory at weak coupling, on a three sphere in the 
grand canonical ensemble with R symmetry chemical potentials. We focus attention 
on near critical values for the chemical potentials, above which the classical theory 
has no ground state. By computing a one loop effective potential for the light degrees 
of freedom in this regime, we show the existence of flat directions of complex dimen- 
sion N, 2N and 3N for one, two and three critical chemical potentials respectively; 
these correspond to one half, one quarter and one-eighth BPS states becoming light 
respectively at the critical values. At small finite temperature we show that the 
chemical potentials can be continued beyond their classical limiting values to yield 
a deconfined metastable phase with lifetime diverging in the large N limit. Our 
low temperaure analysis complements the high temperature metastability found by 
Yamada and Yaffe. The resulting phase diagram at weak coupling bears a striking 
resemblance to the strong coupling phase diagram for charged AdS black holes. Our 
analysis also reveals subtle qualitative differences between the two regimes. 



1. Introduction 



Among the many deep insights to emerge from the AdS/CFT correspondence [T], 
0, 0, is the remarkable connection between large N Yang-Mills thermodynamics 
at strong coupling and black holes in AdS spacetimes. A dramatic consequence of 
this equivalence is that the first order Hawking-Page transition in string theory on 
asymptotically AdS§ x S 5 geometries corresponds to a deconfinement phase transition 
in the Af = 4 supersymmetric, large N Yang-Mills theory on the conformal boundary 
of the spacetime, namely on S 3 x S 1 [|J. Although the above picture most naturally 
applies at strong 't Hooft coupling, it is now well appreciated that such a first order 
deconfinement transition occurs even in the large N free Yang-Mills theory on S 3 x S 1 , 
with thermal boundary conditions |j, El U\. It is believed that the behaviour might 
possibly extend to the weakly interacting theory. 

This raises the exciting possibility that thermodynamics of the weakly coupled 
M = 4 theory on S 3 could be qualitatively similar to the strongly interacting case 
and may provide a window into the physics of black holes in string theory. The 
outstanding issue is to understand how the regimes of weak and strong 't Hooft 
coupling are mapped into one another as the coupling is changed. Of course, it is still 
plausible that, away from the extreme limits of infinite and zero 't Hooft couplings, 
the phase structure at non-zero weak coupling could be qualitatively different from 
that at finite strong coupling. This latter possibility, if true, would imply non- 
analyticities in the theory as a function of the 't Hooft coupling, and qualitatively 
different thermodynamics of string theory on strongly curved backgrounds || [| 
compared to semiclassical gravity. 



For these reasons, mapping out the thermodynamic phase structure of Af = 4 



theory has received much attention from various viewpoints in recent years [|10| , |TT 

ma, [m m m m p, m m S EDJ. 



In this paper we study the weakly coupled Af = 4 theory with SU(N) gauge 
group, on an S 3 of radius R, in the grand canonical ensemble with chemical potentials 
for global R-symmetry charges. The Af = 4 theory has three global U(l) symmetries 
generated by the Cartan elements of the SU(4) R-symmetry. To pass to the grand 
canonical ensemble, we introduce chemical potentials (/ii, \xi, ^3) for the three global 
R-symmetry charges (Ji, J 2 , J3) in the theory. 

We will explore a particularly interesting corner of the phase diagram which 
turns out to be analytically tractable in the presence of chemical potentials. When 
the M = 4 theory is formulated on the three sphere of radius R, the six scalar fields 



- 1 - 



obtain a mass 1/R due to conformal coupling to the curvature of the sphere. In 
addition, the introduction of a chemical potential for a global charge induces an 
effective negative mass squared for scalars carrying that charge. There is thus a 
critical regime of values for the chemical potentials \x v = R~ x + O(X), for which the 
scalars become light degrees of freedom. (Here A -C 1 is the 't Hooft coupling.) We 
will refer to this region as the "near critical" regime. 

The classical theory is only stable when \x v < -R -1 , and this statement is true 
also at one loop order. When the chemical potential exceeds this value, the classical 
Hamiltonian becomes unbounded from below. Note that in a simpler field theory, 
such as a massive complex scalar field theory with a quartic interaction, introducing 
a chemical potential exceeding the mass leads to Bose-Einstein condensation. This 
does not appear to be the case, at least perturbatively, in the Af = 4 theory where 
the quantum corrections at weak coupling are always systematically smaller than the 
tree level mass term. The only situation where there is a possibility for tree level 
and one loop radiative corrections to compete is when fi p ~ R^ 1 + O(X), so that the 
classical potential already contains a term of order A. 

In addition to requiring near-criticality of chemical potentials, for the most part 
we will focus on low temperatures TR <C 1. Our analysis complements the work of 
Yamada and Yaffe who investigated the weakly interacting theory 1 in the tem- 



perature range 1 <C TR < 1/ y/X. One of their interesting results was to demonstrate 
the existence of a metastable deconfined plasma phase, at high temperatures in the 
range R- 1 < \i v < \/XT 2 + R~ 2 . 

In the near critical regime where two or more scalars are light, there exist clas- 
sically (almost) flat directions. They are parametrized by mutually commuting, 
constant background values for the light scalars. At a small finite temperature, the 
shallow directions include mutually commuting constant configurations for the ther- 
mal Wilson-Polyakov line. Moving along these directions generically Higgses the 
theory: SU(N) — > [/(l)^ -1 . For sufficiently large values of the diagonal modes of 
the scalars, the off-diagonal modes of all the Kaluza-Klein harmonics acquire large 
masses. These can be integrated out to generate a Wilsonian effective potential at one 
loop for the light, diagonal modes in the regime of near critical chemical potentials. 

Our first result in the near critical regime at T = is that the one loop radiative 
correction to the scalar potential vanishes, following a non-trivial Bose- Fermi cancel- 



1 The authors of [|TT] also studied the free theory in the grand canonical ensemble, obtaining a 
line of first order Hagedorn/deconfinement transitions for generic values of the chemical potential. 
We do not address the issue as to whether or not the first order line gets modified at non-zero weak 
coupling. 
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lation in regularized Casimir sums in the presence of expectation values for the light 
scalar modes. This means that the nearly flat directions (at T = 0) are not lifted by 
quantum corrections in the vicinity of critical chemical potentials. 

In particular, when the chemical potentials are at their critical values, and at 
T = 0, due to vanishing quantum corrections there is a Coulomb branch moduli space 
of complex dimension N, 2N or 3N, depending on whether we have one, two or three 
critical chemical potentials switched on. For each case we also have a different number 
of zero energy modes: With fii = R _1 ,fi 2 — A*3 = a single holomorphic adjoint 
scalar mode with zero energy appears, while for — iii — -R _1 ,A i 3 = 0, there are 2 
holomorphic scalar zero modes. The situation with three critical chemical potentials 
reveals two adjoint fermion zero energy modes along with 3 holomorphic scalar zero 
modes. 



The appearance of the moduli spaces at critical chemical potential and the asso- 
ciated zero energy modes can be understood more generally as follows. The generator 
of time translations of the M = 4 theory with chemical potentials on S 3 , may be 
expressed as 

3 

A(/x p ) = A-^/i p J p , (1.1) 
p=i 

where A is identified with the dilatation operator of the theory formulated on M 4 . 
With one critical chemical potential /jli = (and fi 2 = A*3 = 0), this operator 
vanishes on all states with RA = J 1; which are the infinite set of | BPS operators in 
the theory. At the critical values for /i p , the superconformal algebra ensures positive 
definiteness of the above Hamiltonian and the | BPS operators thus constitute the 
infinitely degenerate set of ground states of the theory. This can be interpreted as the 
origin of the flat directions at critical chemical potential. For two and three critical 
chemical potentials, the ground states are parametrized by | th and | th BPS operators 
respectively. The dimensions of the moduli spaces we find and associated zero modes 
are consistent with expectations based on our knowledge of the generators of | BPS, 
i th BPS and | th BPS states in the M = 4 theory [0, |23[ . 



At the critical values for the chemical potentials, upon switching on a small non- 
zero temperature TR -C 1 the theory acquires another set of zero modes in addition 
to the light scalars. These new zero modes are the diagonal elements of the Polyakov 
loop matrix. We find a joint effective potential for all the light modes and deduce 
that eigenvalues of the Polyakov loop matrix experience purely a mutual attractive 
force causing them to all collapse on to a point. This corresponds to a deconfined 
phase wherein the trace of the Polyakov loop has non-zero expecatation value. Fur- 
thermore, the theory develops a mass gap due to exponentially small thermal masses 
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at low temperature. This means that thermal effective potential for the scalars has 
a positive curvature near the origin at critical chemical potential. Away from the 
origin, at large field amplitudes it asymptotes to a constant (3/16R). The small 
positve curvature near the origin allows us to raise the chemical potentials beyond 
their classical limiting value 1/R, and still obtain a locally stable configuration at 
the origin. Raising the chemical potential(s) above 1/R, however, causes the scalar 
potential to have a runaway behaviour (unbounded from below) at large field ampli- 
tudes. This results in a metastable state with an exponentially diverging lifetime in 
the N — > oo limit. We expect this low temperature metastable phase to be a smooth 



continuation of the high temperature metastable plasma discovered in [11 



Technically, there is an important difference between the high and low temper- 
ature regimes. At high temperatures, TR ^> 1, the effective potential is obtained 
basically by a flat space computation on IR 3 . The low temperature effective potential 
on the other hand, in the presence of non zero R-charge densities 7^ 0) depends 
on the details of the compact space on which the theory is formulated. Our analyisis 



complements the work of |TT| , filling in the low temperature regime of the phase dia- 
gram of the theory with chemical potentials. The final weak coupling phase diagram 
is shown in Figure 3. 

Perhaps surprisingly, a quick comparison of Figure 3 with the strong coupling 
phase diagram in Figure 4 reveals striking similarities. At strong coupling, the phase 
diagram is dictated by the thermodynamics and stability properties of R-charged 
black holes in AdS space. More detailed aspects of these are discussed in Section 4. 
Here we further remark that it has also recently been found that the region in the 



strong coupling phase diagram, below the black hole instability line and above the 
critical value of the chemical potential, exhibits a metastability . This metastability 
corresponds to a singe (probe) D3-brane splitting from a cluster of large N rotating 
branes whose near horizon geometry is the charged AdS black hole background. 2 This 
is exactly the physics expected from the weak coupling analysis where the lifetime 
of the metastable state is determined by the probability for one scalar eigenvalue to 
tunnel out or be thermally activated into the unstable region. 

Our weak coupling analysis, however, also reveals certain important differences 
with the strongly coupled regime. We find that, in the [i — T plane, the metastable 
region shrinks to zero size at zero temperature - the instability line meets the first 
order deconfinement line at fi = 1/R and T = for any choice of chemical potentials. 
At strong coupling, the black hole instability line and the first order Hawking-Page 
lines meet at \i = 1/R and non-zero temperature. Only for the case with equal 



2 As emphasized in p4j this phenomenon is distinct from the source of local thermodynamic 
instability found in [E5l Era] . 
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chemical potentials do the two lines meet at T = and fi = 1/R at strong coupling. 

The organization of this paper is as follows. In Section 2, we review how R- 
symmetry chemical potentials are introduced in the M = 4 theory. In section 3, we 
show how to compute the one loop effective potential for the light degrees of freedom 
near critical values for chemical potentials. We perform the calculations with one, two 
and three critical chemical potentials at zero temperature. In Sections 3.2 and 3.3 we 
establish the existence of flat directions at zero temperature and their interpretation 
in terms of BPS states becoming light. Sections 3.4 - 3.5 are devoted to establishing 
the existence of the metastable plasma phase at low temperatures. In Section 4, we 
describe the similarities and differences between the phase diagrams at weak and 
strong coupling. Conclusions and future directions are summarized in Section 5, 
and finally, an Appendix is devoted to the spherical harmonic decomposition of the 
theory on S 3 . 

2. R Symmetry Chemical Potentials 

The M = 4 theory has an SU (4) r global symmetry. There are thus three chemi- 
cal potentials which can be introduced associated to the maximal abelian subgroup 
U(l) 3 C SU(4:) R . The six adjoint scalars {(pi} (i — 1,2, ... ,6) transform as the an- 
tisymmetric 6 of SU(4)r, while the fermions are in the fundamental representation, 
the 4 of SU(4:)r. (We follow the conventions of |l]J below). 

We choose the three U(l) generators of the Cartan subalgebra to be 



in the fundamental representation. Packaging the six real scalars into three complex 
combinations, 




\ diag(l, 1,-1,-1), 
|diag(l, -1, 1, -1), 
|diag(l, -1, -1, 1), 



(2.1) 




(2.2) 



we can define the six-vector 



$ = $ 2 , $2, $3,$*). 

The three U(l) generators acting in this representation are then 



(2.3) 



J Rf = idiag(l, -1,0, 0,0,0), 
i$ = ±diag(0, 0,1, -1,0,0), 
i*f = idiag(0, 0,0, 0,1,-1), 



(2.4) 
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which clearly assigns opposite charges to the fields $j and their complex conjugates. 



The grand canonical partition function is defined to be 

Z(T,fi p ) = Tre-^-Zr^) , (2.5) 

where the J p are the three associated conserved R charges. The charge densities 
involve both fermionic and bosonic contributions. With the above choice of the 
U(l) generators, the chemical potential assignments are determined via the following 
expressions 

3 

22 fM R i = iMh — t*2, — A*2, A*3, -A*3) (2.6) 

1=1 

for the six scalars in Eq. (|2.3| ). Similarly, the fermion chemical potentials are deter- 
mined as 

3 

^2 fH R t = diag(/ii, /i 2 , /2 3 , pi) (2.7) 

i=l 

where 



/"I = |(//1 + /i2 + ^3) 
— n,n — Uq") 

(2.8) 



1^2 = |()"1 - A*2 - A*3) 



/i 3 = |(— A*i - A*2 + /i 3 ) 

/2 4 = |(-A*1 + ^2- fa) • 



The grand canonical partition function Eq. ( |2.5|) can also be realized as a Eu- 
clidean functional integral for the theory on S 3 x S 1 . In the functional integral 
or Lagrangean formulation, the chemical potential for each global charge can be 
thought of as introducing by hand, a constant (imaginary) background for the time 
component of a gauge field associated to the respective global U(l) symmetry. In the 
presence of the R charge chemical potentials the Lagrangean for the M = 4 theory 
on S 3 x S 1 becomes 



(3 3 
\F, V F^ + § J](A^2 P -i - ViAo02 P ) 2 + 1 
p=l p=l 



4>2p + i^p5vfi<f)2p~l) 2 + 



2R 2 ^ a V 



(2.9) 



where all derivatives are gauge covariant derivatives on S 3 . The ip A are four compo- 
nent Majorana fermions 

V = M, X Aa f (2.10) 
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The indices a, b = 1, ... 6, and A, B = 1, ... 4, while p — 1,2, 3. The 4x4 
matrices, a p and (3 P are the Clebsch-Gordan coefficients, satisfying 

{a p , a q } = -25 pq , {(3 P , (3 q } = -25 pq , [a p , [3 q ] = 0. (2.11) 

Explicit representations for the a and /3 matrices are given in terms of Pauli matrices 

as, 

«-ct)' «-(-Vo) <-> 

^{~to- ^=u?). ^ 

Classical (in) stability: 

It is immediately clear from the scalar kinetic terms in (|2.9| ) that in the presence 
of chemical potentials the Euclidean action is not real. Notice also that with the 
chemical potential, the conformal scalars have an effective mass given by 

m 2 p = R- 2 -^ p p = 1,2,3. (2.13) 

Hence only when m 2 > 0, does the classical theory have a stable vacuum. For 
jj p > there is a classical instability in the theory along directions in field space 
for which 4>2 P -i and (f>2 P commute. It follows that in flat space, i.e. on M 4 , where all 
the fields are exactly massless, it is not possible to introduce a chemical potential in 
M = 4 SYM since the resulting theory has no ground state (at least classically) and 
the grand canonical ensemble is ill-defined. In finite volume and in particular on S 3 , 
the conformal coupling of the scalars to the background curvature allows for a mass 
which leads to a stable vacuum for a range of values of the chemical potential. 

It is worth noting that this kind of instability in a massive interacting scalar field 
theory, driven by a chemical potential exceeding the mass, generally leads to Bose- 
Einstein condensation, i. e. a VEV for the scalar fields and a spontaneous breaking of 
the U(l) symmetry. For example, this occurs in the (ft 4 theory, where the interactions 
stabilize the ground state at some finite VEV. In the Af = 4 theory, classically at 
least there is nothing to stabilize the theory along the mutually commuting directions 
of configuration space when one or more of the \x v exceed Br 1 . Physically, when 
fi p > B^ 1 the system can always lower the energy by populating the vacuum with 
any number of charged quanta which it can borrow from the bath. 

In the vicinity of the critical chemical potentials fi p ~ -R -1 , the scalar fields of 
the theory are the light, almost massless degrees of freedom. The appearance of these 
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new light modes makes the approach to the critical chemical potential an interesting 
regime to study. In this paper we consider these near-critical regions in more detail. 
There are three cases which will each be considered separately: 

(i) nx ~ BT 1 , ^2 = ^3 = 0; 

(ii) fi^^-B' 1 , 3 = O; (2.14) 



in 



We will consider the approach to the critical chemical, with vanishing as well as small 
non-zero temperatures. 



3. One Loop Effective Potential - Generalities 



In this section of the paper, we set up the calculation of a one-loop effective potential 
for the M = 4 theory on S 3 with a radius R, at both zero and non-zero (but small) 
temperatures in each of the near-critical regions (|2.14j) above. 



Adopting a Wilsonian approach, we will compute the effective potential for the 
lightest degrees of freedom in theory, by integrating out all the other heavy modes 
in the background of the light modes. The natural mass scales for most of the heavy 
degrees of freedom are the inverse radius of the S 3 , namely R~ l and the temperature 
T = All the fields can be expanded in terms of spherical harmonics on S 3 

and Matsubara modes on the thermal circle. The light modes must necessarily be 
constant modes on S 3 x S 1 . The only fields which have such constant modes are the 
scalars (fi a and A , the gauge field component around S 1 . 3 The scalars generally have 
a mass of order since they are conformally coupled to the background curvature. 
However, as we have seen above, in the presence of a near-critical chemical potential, 
[Mp ~ BT 1 , the effective masses of 4>2 P -i and 4>2 P are small and so we should include 
their constant modes in the effective potential. 

Thus there are two sets of light modes in the near critical theory. First, we have 
for each near-critical fi p , light adjoint scalar modes which are the homogeneous parts 
of the respective fields on S 3 : 

<p»-i = 2^3- / 53x51 K-i . = 2^r 3 / 53x5i K ; ih - R- 1 - (3-i) 



3 The other gauge field components A { are vector-valued on S 3 and so do not have a constant 
mode. 
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In addition to these, at finite temperature, we must also account for the spatial zero 
mode of the holonomy of the time component of the gauge field around the thermal 
circle: 

« = f A)- ( 3 - 2 ) 



2n 2 R 3 



The Wilsonian effective potential has a tree-level contribution, as well as loop 
corrections, 

V cS = \- 1 V + V 1 + O(\). (3.3) 

Here A = g 2 N is the 't Hooft coupling. Although we only restrict to one loop 
computations, we will consistently express all quantities as fucntions of the 't Hooft 
coupling and N, indpendently, since these are the parameters relevant for evenetual 
comparison with the large N gravity dual. The tree- level term is 

H = ^^Tr(-^[A o ,0j 2 -^[0 a ,0 b ] 2 + ^(i?- 2 -^)(0^_ 1 + 0y) . (3.4) 

a a<b p 

When the chemical potentials are small so that fi p < -R -1 , the tree-level potential 
forces the associated scalar VEVs to vanish. On the other hand, if fi p > R~ l , 
the theory becomes unstable. In the region of one or more near-critical chemical 
potentials, fi p ~ -R -1 , the classical potential has almost flat directions. These flat 
directions correspond to mutually commuting VEVs f2 P -i and v?2 P , in addition to a 
mutually commuting a. 



One loop correction near-critical chemical potentials 

We will now outline the computation of the one-loop effective potential in the theory 
approaching the critical chemical potential in each of the three regions ( |2.14| ). The 
near-critical region for the chemical potentials will be defined specifically as 

li p — R~ l + 0{\). (3.5) 

This parametrically small approach towards criticality is chosen so that the tree-level 
potential for the commuting modes is of the same order as the first loop correction. 
It is in this situation that there is the possibility of competition between these two 
effects with potentially interesting physics. If — R~ 1 \ is parametrically larger, 
then depending on the sign, the theory will simply either have a stable vacuum at 
the origin, or have perturbatively unstable runaway directions. 

For each critical chemical potential /i p , classically the theory will have almost flat 
directions parametrized by (a holomorphic combination of) the diagonal elements 
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of mutually commuting matrices ip 2p -i and (p 2p , (p = 1,2,3). In addition, there 
is another set of moduli parametrized by the diagonal elements of a, which also 
commutes with (p 2p and (p 2p -i at a minimum of the classical potential. Therefore at 
a generic point along this flat potential we have, 

a = diag(ai) <p 2p -i = diag (</?2 P -ii) , ip 2p = di&g(ip 2pi ) , i = 1, 2, . . . , N. 

(3.6) 

Recall from ( ft.ip that the fields (a, (p a ) are the background values for the spatially 
homogeneous parts of the full quantum fields (Aq, <j) a ). 

Near-critical chemical potential for only one of the three U(l)s 

We begin our analysis with the simplest situation where only one near critical chem- 
ical potential \i\ is turned on: 

fit ~ R" 1 ,^ = fJ-3 = 0. (3.7) 

In this situation, the diagonal elements of the fields <pi,<p 2 and a are the lightest 
modes in the theory. 

At generic points of the classically flat directions parametrized by the diagonal 
elements of (pi and (p 2 , we will integrate out all the inhomogeneous modes (the 
Kaluza-Klein harmonics on S 3 ) as well as the off-diagonal homogeneous fluctuations. 
In order to ensure the validity of the semiclassical or one loop approximation, it will 
be necessary that the off-diagonal fluctuations have relatively large masses. From 
the Lagrangian ( |2.9| ), we can estimate the masses of the light off-diagonal excitations 
of the zero momentum scalar modes to be 

™% ~ + #y + R- 2 - fi\ (3.8) 

where tp a>i j = ip ai — <p a j. Validity of perturbation theory in the Wilsonian sense, 
about diagonal scalar backgrounds then requires that the off diagonal modes be 
much heavier than the light diagonal degrees of freedom: 

R2 E (via + via) » I 1 - (^) 2 I ^ 0{\). (3.9) 

p 

This condition is easy to ensure in the near critical region ( |3.5| ) by choosing appro- 
priately large VEVs for all the diagonal entries, such that the differences between 
them also remain parametrically large. 

The masses of all heavy excitations will depend only on the sum 

<Pij= (flij + ^hj) , (3-10) 
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as seen above. 4 For this reason, at leading order in the coupling, the one-loop contri- 
bution to the effective potential for light modes only depends on this combination. 
In this case, we can use the resulting symmetry to only turn on a VEV for one of 
the scalar fields, which we can take to be pi, for instance. The full dependence on 
the scalar fields can then be reconstructed by replacing </??• (we drop the subscript 



The radiative corrections at the one loop level are then obtained by taking the 
constant modes as backgound VEVs and integrating out all the massive modes of 
the fields. To this end, we shift (pi — ► (pi + p and at non-zero temperature, we also 
introduce a zero mode for the time component of the gauge field via Aq — > Aq+oi. The 
VEVs contribute to the effective masses for the modes and the one-loop correction 
involves the logarithm of the resulting fluctuation determinants and hence depends 
on the VEVs in a non-trivial way. It is important to note to leading order in the 
coupling in the near critical region, we can take fii = R^ 1 exactly. 

Before we compute these determinants, we must first fix the gauge. We find 
it convenient to do this by working in a conventional gauge of a spontaneously 
broken gauge theory and then to specialize to Feynman gauge £ = l. 5 We add to 
the action the gauge fixing term 



In the above equation and in what follows, we leave adjoint action by <p as implicit, 
i.e. pep = [<p, <p], p 2 (p = [ip, [<p, 0]], etc. In addition, D = <% + ia includes a the zero 
mode part of A only, and, as for ip, adjoint action for a is implied. 

Although in the absence of chemical potentials the gauge fixing removes cross 
terms between the gauge field and the scalars, the presence of a chemical potential 
and a VEV <p introduces additional ones that are not removed by the gauge fixing 
term. The modes A , <pi and (p2 are all coupled together as seen from the expression 



4 In the presence of diagonal background VEVs for the lightest scalars, all the degrees of freedom 
which are integrated out are necessarily the off-diagonal fluctuations of every KK harmonic on S 3 . 

5 We shall not be unduly concerned by the usual £ dependence of the effective potential that 
generally plagues gauge theories. The reason is that in the vicinity of the critical region there is an 
exactly flat direction at tree level and so the VEV can be taken to be non-zero whilst remaining on 
shell. The £ dependence then drops out as one can explicitly find by including it in all subsequent 
steps. 



1) by O. 




(3.11) 
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for the bosonic part of the action at Gaussian order in quantum fluctuations: 



^(boson) = 1 Tr h Ao (-D 2 - A« + V 2 )A + \M-Dl - + v 2 )A t 



9 2 



2 ^uv —0 — 1 T /--HI i 2 

+ |0i(-,Do 2 - A (S) + V 2 )0i + |02(-^ o 2 " A (S) + V 2 )02 



6 

2\ 



+ c(-D 2 - A« + <p 2 )c 



a=3 



(3.12) 



Here we have set /ii = -R -1 , A^ s ^ and are the scalar and vector Laplacians on 
S 3 , and we have explicitly included the ghosts c, c in the Lagrangean. The detailed 
properties of these Laplacians and their eigenvalues are summarized in Appendix A. 
We need only to note that while A* 5 * 1 has an £ = mode, the vector Laplacian A^ 
on S 3 does not have a zero mode. Furthermore, the vector fluctuations Ai can be 
decomposed into the image and kernel of the gradient operator V* as A 1 = B l + C l 
with with V i B i = and G i = V 1 f . 

The fermionic fluctuations about the non-zero backgrounds for (pi and if2 are 
governed by the Lagrangean, 



^(fermion) = Tr ^ j Q _ ^ _ ^ + j^^) ^) , ( 3 . 13 ) 



where 



Ha — R (|,|)— !' — D' (3-14) 



using ( p.8|) . The ^ are four component Majorana fermions. 

The one loop contribution to the effective potential is obtained by integrating 
out each massive fluctuation, giving rise to the corresponding fluctuation determi- 
nant. Importantly, the quadratic fluctuation operators for the (A , 0i, 2 ) sector and 
the fermionic sector are off-diagonal in flavour space. For example, using that the 
eigenvalues of A^ s ^ are given by £(£ + 2)R~ 2 , for £ = 0, 1, . . ., we can evaluate the 
determinant of the fluctuation operator coupling the (A , $1, </> 2 ) sector for a given 
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spherical harmonic number, 

(-D 2 + £(£ + 2)R~ 2 + y? 2 2iTV 

do! -D\ + £(£ + 2)R~ 2 + if 2 2iR~ 1 b 

-2i?'V -2iR- l D Q -D 2 + £(£ + 2)R- 2 + ip 2 J 

(-D 2 + £(£ + 2)R- 2 + V 2 )(-D 2 + £ 2 iT 2 + <p 2 )(-D 2 + (£ + 2) 2 R- 2 + y? 2 )] . 

(3.15) 

The zeroes of this expression viewed as a polynomial in D yield precisely the energies 
of the harmonics. Note that for this, the simplest situation at hand, these match the 
results summarized in Table 1, if we set ji\ = , /i2 — A*3 = 0. The table shows 
that introducing a chemical potential leads to a re-organization of the energy levels. 
In the bosonic sector, when ^ = R~ l is turned on, the result is three new towers of 
modes that are shifted in such a way that one of the new towers is identical to the 
original A tower (which can be cancelled with the Cj and ghosts) while the other 
two new towers are the deformations of the original <p\ and 02 towers. For fermions, 
notice that the original half-integer graded modes are now integer graded. 

The fluctuation determinant for each species then enters the effective potential 
at one loop as 

T i N °° 

V ^ = ^5^2 £ logdet[-D 2 + ^-) 2 ] , (3.16) 

species ij=l t=£o 

where df^ is the degeneracy of bosonic (fermionic) modes with angular momentum 
quantum number £. The integer £ Q is the lower limit on the angular momentum 
quantum number. The quantity ee is the energy of the mode in question. The 
degeneracy factors d e are positive or negative depending upon the the statistics 
of the corresponding fields. 

In Table 1 we summarize the data associated to each set of modes. We emphasize 
that these energies are strictly only correct when either /j, p = or fi p = it! -1 . Where 
sign choices exist all the possible combinations must be taken. 

The contributions to the effective action are standard expressions in thermal field 
theory. First of all, the eigenvalues of id are 27rn//3, for n e Z for bosons. When 
acting on fermions, due to antiperiodic boundary conditions for fermions around the 
thermal S 1 , the operator id has eigenvalues 2n(n + 1/2)//?, for n G Z. It is standard 
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1 
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C 


(7 4- ll 2 
I* T 


,/B- 2 £(£ + 21 4- 09 2 

-wjl t^-c T T ^ 


1 


(c, c) 


-2(£ + l) 2 


Y^" 2 ^ + 2) + ^ 2 





(A),01,2)l 


(£ + l) 2 


l/R~H{l + 2) + y>2 





(^0,01,2)2,3 


(£+l) 2 


v / i?- 2 (^+l±i?/ii) 2 + ^ 2 





03,4 


(£+l) 2 


v / i?- 2 (^ + l) 2 + ^ 2 ±/i 2 





05,6 


(£ + l) 2 


Y / i?- 2 (^ + l) 2 + ^ 2 ±/i 3 







-£(£+1) 


^R~ 2 {£ + \ ± \R^i) 2 + <^ 2 ± f ± f 


1 



Table 1: The fields, their degeneracy and energies as a function of the chemical potentials 
with a nonvanishing VEV for (p\. The expressions for the mode energies are only valid 
when either fj, p = or \i v = i? . 

practice in thermal field theory to perform a Poisson resummation over n in such 
a way that each contribution splits into a piece that describes the theory at T = 
and the non-trivial "thermal" part which vanishes as TR — > 0. The zero temperatue 
piece is the Casimir energy in the presence of background expectation values for 
fields. So a typical term in the one loop potential (|3.16|) can be expressed as 



. 00 

— — — ^rf £ log det(-D 2 + e e (cp) 2 ) 



Bosonic : 



. N OC _. 00 . 

y^) 2 E E «<f " i E ie-*(«)l 

V ' ij=li=i H n=l 



N oo , oc , -. 



1 1 r 1 ( — 1 ) n "i 

Fermionic : vous^)2 - ^E^ e_ri/3N( ^ )lcos (^)} 

(3.17) 

where we have written the adjoint trace explicitly and denned <pij = tpi — tpj and 



More than one non-zero critical chemical potential 

The analysis with more than one non-zero critical chemical potential proceeds sim- 
ilarly to the above. The only difference lies in the mixing matrices for bosonic and 
fermionic fluctuations and their resulting eigenvalues. 
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When [Mi ~ 1^2 — R' 1 and [13 = 0, the light scalar modes are (<pi, . . . ^4) and a 
5x5 mixing matrix results for the fluctuations in the bosonic sector. The zeroes of 
the determinant of this mixing matrix in frequency space correspond to the mode 
energies. Similarly one may find the fluctuation energies for the case of three critical 
chemical potentials when all six scalar fluctuations and Aq are coupled. The results 
for all the mode energies, both bosonic and fermionic are listed in Table 1. 



4. The zero temperature effective potential 



As we take the low temperature limit of ( |3 . 1 7| ) , the thermal contributions are expo- 



nentially suppressed by the Boltzmann factors. Then we are only left with a Casimir 
energy sum that yields the one loop effective potential 

N 00 

Vi{T = 0) = voi(M2 £ (4-1) 

Species ij=l i=ia 

As is well known, this sum is formally divergent and needs to be regulated. 

The choice of regulator is a subtle issue especially in the presence of chemical 
potentials and there is more than one way to regulate such sums. We choose to cut 
off the sums keeping two crucial points in mind: First, the cutoff will be imposed on 
the energies of the modes rather than the angular momenta i so that the regulator 
is general coordinate invariant. 6 Secondly, the energy cutoff function will be chosen 
to be completely independent of the chemical potential(s). What this means is 
that, all mode sums will be regulated using the cutoff functions of the theory with 
fix = Hi = // 3 = 0. This requirement may be motivated by the physical observation 
that introducing a chemical potential is a deformation of the state of the theory and 
not the functional integral measure. 

For each field type, we introduce the regulated expression for the energy sums 

N 00 

4* =2^E d ? (F) M^)l/(4 0) /A)- (4-2) 
ij=i e=e 

Here f(x) can be thought of as the smooth version of a cutoff function which is 
unity for x < 1 and for x > 1 such that /(0) = 1 while all the derivatives 
f'(0) = /"(0) = • • • = 0. Crucially, the cutoff is on the energy ef^ which is defined 
to be the energy of the respective mode with \i p = 0. 



3 See the footnote (30) in H and Chapter 6 of [p7j for a discussion of these issues. 
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We can evaluate the regulated sum by using the Abel-Plana formula [28] appro- 



priate to a function with branch cuts on the imaginary axis or the left half plane 

E F H=/ dxF(x) + -F(0)-2 dx-^-. (4.3) 
n=0 J ° Z Jo e -1 



Let us now turn to the evaluation of each field contribution to the Casimir energy 
seriatim. 

The first point that is immediately clear is that large cancellations occur between 
the Ci, the ghosts (c, c) and (A o ,0i,02)i which is basically the A field. Indeed, all 
the £ > contributions cancel between these, leaving a net result — l\(f\- 

The regularized contributions from the Bi and the scalar fields remain in the 
bosonic sector. Taking care to use the \i v = energies in the cutoff functions and 
applying the Abel-Plana formula, we find that all bosonic modes together induce the 
one loop potential 

S B =A 4 R 3 - -RA 2 - R 3 ^ 2 A 2 + — - ~<p*R 
2 v 12R V 



1 (WW A \ r x 2 Jx 2 R~ 2 -if 2 

+ -^R 3 log + c ' V 



(4.4) 



2 r ~° V 2A J 1 " J Rv e 2 ™ 



Interestingly, the potentially problematic linear term — ||v?| from the ghosts, A$ and 
the Ci, cancels off in the full sum. The result is valid for any number of critical chem- 
ical potentials. The sum over energies for more than one critical chemical potential 
reduces to the one for a single chemical potential due to a simple cancellation of the 
dependence on /i 2 and This cancellation is obvious from the fluctuation energies 
of 03, ... 06 listed in Table 1. 

For the fermions, again employing a cutoff on the mode energies at zero chemical 
potential, the Casimir energy is 

$f = -a 4 r 3 + Ira 2 + rVa 2 + 4b + -^R 

2 48R 4 

(4.5) 



2 V 2A J J Rip e 2 ™-l 

The total one loop quantum effective potential on S 3 is therefore 

' (*F + *b) = t4f7^ 4b- ( 4 -6) 



Vol(5 3 ) v r Vol(5 3 ) 16R' 
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This is remarkable since all the dependence on the background VEVs has completely 
cancelled out to yield the zero point energy of M = 4 theory on S 3 . This implies 
that at least at one loop order, at the critical chemical potential, there is a complete 
cancellation between bosonic and fermionic quantum fluctuations resulting in flat 
directions in the quantum effective potential. 7 

To summarize the result of the computation presented in this section: We have 
found that at the critical value for one, two or all three chemical potentials of M = 4 
theory on S 3 at zero temperature, there is no radiatively induced potential for the 
scalar fields at one loop order. This implies, at least at the one loop order, a Coulomb 
branch "moduli space" of complex dimension N, 2N or 3N, depending on whether 
one, two or three chemical potentials are at the critical value. 

Below we offer a natural explanation for the appearance of these flat directions 
in the effective potential calculation, and its existence to all orders in the interacting 
theory at zero temperature. 



4.1 Massless BPS states at critical chemical potential 

At zero temperature, chemical potentials for the three U(1)r charges deform the 
Hamiltonian on S 3 as, 

3 

A^A-^/i p J p . (4.7) 
P =i 

The M = 4 superconformal algebra p2| includes the following anticommutation rela- 



tion between the superconformal generators Sa and the supercharges Q B/3 . Noting 
that on S 3 , following the rules of radial quantization, (Q Aa )^ = S^a, we have 

{Q Aa \Q^} = 6i6^A + J2(J P R P f A €+ Si • (4-8) 



A similar relation holds for Q and S. The indices A, B = 1 ... 4 are £77 (4) # indices 
and Rp (p = 1, 2, 3) are the three U(1)r Cartan generators in the fundamental rep- 
resentation of the SU(4:)ii algebra; J\ is a generator of an SU(2) C £0(4) isometry 
group of S" 3 . 

For a subset of the supercharges, the right hand side of the above relation vanishes 
on BPS states of the N = 4 theory. Multi-trace operators with R charge (Ji,0,0) 



7 It is worth contrasting this result with the same calculation done at zero temperature and 
vanishing chemical potentials. In that case one obtains a complicated one loop effective potential 
on S 3 . 
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and with (-RA) — J\ = are \ BPS states which are annihilated by one half of the 
chiral supercharges in (|4.8| ). Similarly, states satisfying (RA) — Ji — J 2 = and 
(RA) — Ji — J 2 — J3 = correspond to j th BPS and | th BPS states respectively. 

The operators A - Ji/R; A - ( J x + J 2 )/R and A - (J x + J 2 + J 3 )/R are of 
course the Hamiltonians of the theory with critical chemical potentials ( (4.7[ ). It is 
thus clear that the ground states of the Hamiltonian with critical chemical potentials 
are infinite sets of one half, one quarter and one-eighth BPS states of the TV = 4 
theory, depending on whether we have one, two or three critical chemical potentials. 
Importantly, in each case there is an infinite degeneracy of ground states, all of which 
have zero energy, following from 



Choosing an appropriate Af = 1 subalgebra of the J\f = 4 superalgebra, the chiral 
ring is the set of holomorphic gauge invariant operators made up of polynomials of 
$1 = 0i + i(p2, $2 = 03 + 204 and $ 3 = 5 + i0 6 , and the chiral gauge field strength 
W a , modulo F-term relations (on M 4 ). These are the | BPS states. On M 4 , the 
chiral primary operators at a generic point on the Coulomb branch of the theory, 
are the the same as the operators of the chiral ring described above. The F-term 
conditions ensure that the fields may be simultaneously diagonalized and the chiral 
operators are distinct polynomials of the 3iV bosonic and 2N fermionic diagonal 
eigenvalues, invariant under the permutation group Sn- On S 3 , these operators are 
built out of the s-waves of the diagonal elements of the three complex scalars and 



two polarizations of the i = 1 harmonic of a fermion [29 . 



Similarly, the half BPS states correspond to Sn invariant polynomials involving 
the diagonal elements of the s-wave of a single bosonic holomorphic field, namely 
<fi + iip 2 . The quarter BPS states are generated by diagonal elements of s- waves of 
the holomorphic bosonic operators <pi + iip 2 and y?3 + i(p±. 

The dimensions of the moduli spaces encountered at the different numbers of 
critical chemical potentials, and the corresponding counting of zero modes is therefore 
perfectly consistent with the interpretation that at critical chemical potential BPS 
states of the M = 4 theory become light. Each point on the resulting moduli space 
represents a coherent superposition of the light BPS states. 



5. Low temperatures and a metastable plasma phase 



To summarize the story thus far, the phase diagram of the free theory was studied 
in P"T| at all temperatures and chemical potentials. At infinite N, the authors of 




x = o 



Figure 1: The A = phase diagram as a function of temperature and chemical potential, 



found by [11 



11] established the existence of a first order Hagedorn/deconfinement transition for 
any fixed \i p as the temperature is increased (see Figure 1.). For n v > R' 1 , the free 
theory ceases to have a well-defined ground state in the grand canonical ensemble. 



In the vicinity of \x v ~ -R -1 , we have shown that light scalars appear and remain 
light in the interacting theory at T = 0. Furthermore, we see flat directions associated 
to these light modes. For a small non-zero temperature TR C 1, we expect the scalar 
modes (<P2p-i,f2p) to continue to be light degrees of freedom in the near critical 
region. In addition there will be another set of zero modes in the theory, namely the 



the diagonal elements of the Polyakov loop matrix. 



5.1 Low temperatures and /Ui = R 1 ,fi2 — ^3 — 



At the critical chemical potential, and at zero temperature, the effective potential 
is flat including quantum effects. There are no non-analyticities in the effective 
potential near the origin, even when off-diagonal modes become light. This is due 
to an exact cancellation between the Vandermonde measure factor (equivalently the 
ghosts and A ) and the contribution from the new scalar zero mode which appears 
at the critical chemical potential. 



Our strategy will be to stay at H\ — R 1 , n% = A*3 = and switch on a small non 
zero temperature TR <C 1. The finite temperature effective potential at this critical 
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point is 



V + V 1 = 

N 1 
Vo\(S 3 ) 



16R 



oo oo " 

fc=l €=1 



The effective potential above is a function of the light fields 



(5.1) 



(5.2) 



A notable feature of this expression is that there are no non-analytic terms near 
the origin where we expect off-diagonal modes of 4>i and 4> 2 to become light. The exact 
cancellation between the scalar zero mode and measure factors persists at non-zero 
temperature. This is related to the fact that with one critical chemical potential, the 
ground states at zero temperature are parametrized by multi-trace gauge-invariant 
combinations of the s-wave of a single holomorphic field, namely tpi +itf2- All other 
modes are heavy. At the origin, therefore, we expect that even if off-diagonal modes 
of the field appear to become light at non-zero temperature, we may gauge rotate 
them away completely. (By the same logic, we expect not to be able to do this for 
the | th and | th BPS states, since they involve multi-trace combinations of more than 
one holomorphic field.) 

The radiative corrections vanish exponentially at large field amplitudes \tpij\R ^> 
1. This is actually a generic, robust feature expected of the effective potential of 
M = 4 theory on S* 3 . In the limit of large field amplitudes \tpij\R ^> 1, all off- 
diagonal modes are proportionately heavier and therefore should decouple from the 
theory. That this happens unambiguously is linked to the UV finiteness of the j\f — 4 
theory. Furthermore, as already noted above, the putative scalar zero energy mode 
(near </??■ ps 0) has cancelled against the ghosts. 

Now we need to determine the minimum of this joint potential. It is clear that 
the Oii experience a purely attractive pairwise potential near tp^ = 0. This is because 
the repulsive Vandermonde determinants have been eliminated through cancellations 
with the zero mode. Hence at finite temperature, at the critical chemical potential, 

o,i — a-2 — . . . — otN'-, (5-3) 

which means that the large theory should be thought of as being deconfined and 
the large distribution of the at is a delta function. With a i3 - = 0, it is clear that 
at <fij = 0, the scalars have a mass squared which is strictly positive and the origin 
is the stable vacuum of the theory. 
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Expanding the potential for R 2 (pfj "C 1, and TR <C 1, we see that thermal effects 
contribute a small mass to the scalars 

V ° + Kl W 2^ ( 16 ~ 8TR e ~™ + 2 ^ 6 ""iV 5>" + ^ J • (5 ' 4) 

Interestingly, a small temperature TR <C 1 has lowered the vacuum energy at the 
origin to 3/16/2 — 8Te -1 / TR and the light scalars <p\ and y?2 have acquired an expo- 
nentially small positive mass. 



Metastable plasma for ji-y > R 1 ,/i 2 = Ha = 



We now want to argue that since a small temperature makes the light scalars massive 
at Hi = R^ 1 , the origin will continue to be a locally stable vacuum if we raise /xi by a 
sufficiently small amount above its critical value. This is of course only possible if the 
thermal mass can beat the instability induced by the chemical potential. Crucially, 
we also require that, for our analysis to apply 

IM X -R- 1 <0(X), (5.5) 

which allows us to set Hi = R" 1 in the one loop calculation, the error in doing so 
appearing at one higher order in perturbation theory. 



With this choice of parameters, we have 



ij = l 



16R 



- 8T e ~T^V 1+R2 v 



(5.6) 



Since the quantum correction vanishes exponentially at large field amplitudes, the 
potential will revert to its classical runaway behaviour in these regions whenever 
/ii > R" 1 . Close to the origin, however, the situation is different. For the following 
range of values of /ii 

0< Ail - J R- 1 <^e-^, (5.7) 

the theory has a metastable vacuum at the origin. The behaviour of the effective 
potential is shown in Figure 2. 



It is reasonable to assume that this metastable phase is the continuation to low 
temperatures, of the high temperature metastability found in |l I] . As in the high 
temperature situation, the lifetime of the state is determined by the probability for 
one diagonal mode, say tpu, to make its way out into the unstable region. This can 
happen either by thermal activation over the barrier, or by tunnelling through the 
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Figure 2: The low temperature TR <C 1 effective potential for one diagonal mode of the 
light scalar with < [L\ — R^ 1 < (4X/tt 2 R) exp(— ^) and ^2 = ^3 = 0. 

barrier. The probability for thermal activation can be estimated by computing the 
Boltzmann suppression factor. At low temperatures , the height of the barrier is 
exponentially suppressed and oc exp(— so that the associated probability for 
thermally exciting one eigenvalue goes as exp (— iVe -1 / TjR ). Estimating the barrier 
penetration probability is more complicated due to the form of the effective potential, 
but the dependence on iV is obvious, due to the scaling of the action with N. The 
lifetime of the metastable plasma diverges in the strict limit iV — > 00, as was also 
found for the high temperature metastable phase in [|ll|] . 



5.2 Two and three non-zero critical chemical potentials 

We have seen that with two and three critical the theory at T = has flat 
directions which are not lifted by quantum corrections. 

At finite temperature however, the one-loop effective potential for the diagonal 
modes with two and three critical chemical potentials displays some qualitatively 
different features compared to the case with one chemical potential. The main dif- 
ference is due to the appearance of additional zero modes for small (fij. The potential 
for the diagonal modes exhibits non-analytic behaviour near the origin. We may un- 
derstand these non-analyticities as being due to the off-diagonal modes of the s-wave 
components of the light holomorphic fields getting excited near the origin, at finite 
temperature. 

At zero temperature and critcial the light Coulomb branch states consist 
of | th and | th BPS states. However, close to (p^ = at finite temperature, there 
are also light off-diagonal excitations which are non-BPS. Near the origin, therefore, 
it is more appropriate to look at the perturbative result for the effective potential 
for the s-wave of the full matrix fields obtained by integrating out higher harmonics. 
As we will see below the only effect of this is to provide a small positive thermal mass. 
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Two critical chemical potentials 



The low temperature analysis (TR <C 1) for two near critical chemical potentials 
yields the following one loop correction to the effective potential for the diagonal 
modes, 



1 ^ Vol(,S 3 ) 116R 2 ^ n \ T ) 

11=1 v ' n=l 

(5.8) 

_ 2 Te-^V 1+ ^ R2+ WR C os + 0(e-™)] . 

Here = X]a=i(v 9 aj — Vaj) 2 - Once again, near (/?^ = 0, the ccj experience a purely 
attractive potential so that ojj = a, and the theory is in the deconfined phase. This 
expression for the effective potential exhibits non-analytic behaviour at the origin. 
With otij = 0, we have 



N 

V x = V - \— + - log ( 1 - e-^v/i+^+^^l 

^ Vol(S 3 ) Ll6i? + 2 8 V J 



(5.9) 



— 2T e^^V^j^+^fi + 



Near the origin, the logarithm provides an infinite attractive potential and is the result 
of integrating out off-diagonal elements of the £ = harmonic of (p 3 +i(p 4 . The second 
term which is 0(e^ 1 ^ 2RT ), originates from integrating out the next lightest state, 
namely the fermion with mass 1/2R. Expanding this contribution about 0^ = 0, 

2 Te~^V T +^+2Wr ~ — L_ ^_2T e - ^ + R ^ ? (5.10) 



Vol(,S 3 ) 2ir 2 R 3 

we see a thermal mass for the diagonal modes. By gauge invariance at the origin, 
the off-diagonal fluctuations will also obtain the same thermal mass. The resulting 
picture is therefore rather similar to what we have seen earlier. Instead of integrating 
out light off-diagonal modes near the origin, we may simply compute the thermal 
mass of the lightest fields, namely tpi + iipi and tp^ + iip 4 by computing the usual 
Feynman one-loop self-energy graphs at finite temperature. At quadratic order, the 
thermal effective potential for these fields close to the origin will have the form 

N 2 ( 3 1 1 A 1 , 9x \ 

V = [l6R ~ 2T 6 " W + 2 R ^ iv Tr ^a) + qnarticj . (5.11) 
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The higher order interactions are the terms responsible for inducing non-analyticities 
in the effective potential for the diagonal modes when off-diagonal fluctuations are 
(wrongly) integrated out close to the origin of field space. 

Hence, a finite temperature introduces a small dip in the potential energy at 
the origin, which then asymptotes to the constant value of 3/(16i?) at large field 
amplitudes in accordance with ( |5.8| ). 



Three critical chemical potentials 



The situation with three critical chemical potentials is similar to that with two 
critical chemical potentials. There is an additional light scalar as well as a fermion 
zero mode. The potential for the light diagonal modes is 

^ VVol(S 3 ) L16.R f^2k-l V J TJ 



(5.12) 



As before the oij experience an attractive potential resulting in a delta-function dis- 
tribution at large N for these fields. There is then a logarithmic attractive potential 
between the light scalar diagonal modes. It is therefore more appropriate to keep the 
all s-wave fluctuations of the matrices <fi,<f2, ■ ■ -fe while integrating out all higher 
partial waves on S 3 . This leads to a temperature induced positive curvature in the 
effective potential near the origin 



N 

V 



2ir 2 R 3 \ 16R 



6 \ 

12Te~^ +8e~^J]- Tr(p 2 a + ... J . (5.13) 

a=l ) 



At large field amplitudes the effective potential approaches a constant, given by the 
Casimir energy of the M = 4 theory on the sphere. 



Metastable phase 



Since we have established a positive curvature for the effective potential at the origin, 
the argument for the existence of a metastable plasma phase above critical chemical 
potential proceeds exactly as in the situation with one chemical potential. 
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For two equal chemical potentials, the metastable phase exists provided the larger 
of the two chemical potentials exists in the range 

2A i 

< ti-RT 1 < — - e~^R. (5.14) 



For three chemical potentials, the metastable plasma phase exists if the largest 
of the chemical potentials satisfies 

8A 

< fi-R' 1 < —- e'TK. (5.15) 
n 2 R 

In both cases, the lifetime of this phase is determined by the probability for one diag- 
onal mode to make it out into the unstable region. This probability is exponentially 
suppressed in the large N limit. 




Figure 3: The known features of the weak coupling phase diagram of N = 4 theory with 
a chemical potential. 



5.3 Metastable plasma at high temperatures 



The picture that we have found above at low temperatures and near the critical 
chemical potential, nicely complements the results of [IT|. In [IT|, a careful high 



temperature analysis of the theory with chemical potentials on S 3 was performed. 
The main observations therein maybe summarized as follows. At high temperatures 
TR 3> 1, the theory behaves much like the theory on flat space. One may then derive 
the "electric" effective theory by integrating out all non-static fluctuations and all 
fluctuations on length scales of order 1/T. This theory is valid on length scales of the 
order of the Debye or electric screening scale ~ (\/AT) _1 . Specifically this means that 
the S 3 radius R ~ {\f\T)~ 1 and is parametrically smaller than the nonperturbative 
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magnetic screening scale. The electric effective theory has an effective mass for the 
scalar fields, given by 

m\ = R~ 2 - $ + AT 2 ; p= 1,2,3, (5.16) 

where AT 2 is the thermal or Debye mass for the scalars. This means that the effective 
potential for scalars will have a positive curvature at the origin, if the largest of the 
three chemical potentials satisfies 

Hp < VR~ 2 + AT 2 . (5.17) 

At large field amplitudes, perturbative quantum corrections are expected to vanish 
due to large masses for the heavy states which are integrated out. In this regime 
the effective potential will revert to its classical behaviour. If the largest chemical 
potential exceeds 1/R, then the potential has runaway behviour for large field ampli- 
tudes. Thus we are left with a metastable vacuum at the origin at high temperatures. 
The lifetime of this vacuum diverges exponentially at large iV and as in the cases 
encountered in this paper, it is determined by the probability for one eigenvalue to 
make it out to the unstable region by thermal activation or barrier penetration. At 
high temperatures TR ~ 1/vA this probablity is oc exp(— N(TR) 3 ) ~ exp(— iVA~ 3 / 2 ) 
which vanishes exponentially in the strict large N limit. 



Putting together the high T results of |TIJ and the low T results in this article, 
we obtain a phase diagram for the theory at weak coupling, summarized in Figure 3. 
(It is, of course, not known as of now whether the first order line seen in the A = 
theory, persists in the interacting theory.) As we discuss below, this weak coupling 
phase diagram bears a remarkable resemblance to the one at strong coupling (Figure 
4). However, there also appear to be certain intriguing and unexplained differences 
which we describe in some detail below. 



6. Comparison with strong coupling 



In the strong coupling limit A — > oc, the SU(N), M = 4 theory at infinite N, 
with R symmetry chemical potentials is dual to Type IIB supergravity in AdS& x S* 5 
with three independent angular motions along the internal S 5 directions. Work- 
ing in global coordinates, there are two possible saddle point configurations that 
contribute to the semiclassical gravity partition function. One of these is the "spin- 
ning" AdSs x S 5 geometry while the other is a charged black hole in AdS. In the 
case where all three rotation parameters are equal (corresponding to three equal 
chemical potentials), the gravity dual reduces to Einstein-Maxwell theory on AdS$ 
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linear 



X - oo ^ 

Figure 4: The phase diagram of N = 4 theory with chemical potentials at strong coupling 
and infinite N . The figure corresponds to a generic pattern of chemical potentials, or 
angular velocities along the internal S 5 of the gravity dual. It has been argued |^4| that 
the region above \i = 1/R and below the instability line is actually metastable in much the 
same way as the field theory at weak coupling. 



space. The system undergoes a first order Hawking-Page transition, as the temper- 
ature is increased for fixed chemical potential, from spinning thermal AdS space to 
the Reissner- Nordstrom- AdS black hole A similar phenomenon occurs for more 
general rotation parameters wherein the dual geometry is described by charged black 
hole solutions in 5D M = 8 supergravity JJl], ^6 . 



As is apparent from the strong coupling phase diagram in Figure 4, fi — R" 1 is 
not a point of instability, except possibly at T = 0. Indeed, at T = 0, setting the 
chemical potential equal to corresponds to the internal S 5 rotating at the speed 
of light. Crossing this limit triggers an instability due to the time direction becoming 
spacelike near the origin of the space . 



When the theory is in the charged black hole phase (deconfined plasma), in- 
creasing the chemical potential beyond 1/R does not lead to any local instabilities. 
Local instabilities are only triggered at much higher values of the chemical potential 
The fate of the theory above this instability line remains unknown. 



Figures 3 and 4 illustrate the remarkable qualitative similarities between the 
weakly coupled (one loop) field theory on the one hand and its strongly coupled 
limit on the other. In fact it has recently been argued |23] that the metastable 



plasma phase seen at weak coupling has a strong coupling analogue and that the 
charged AdS black holes are metastable for fi > R~ Y and below the instability line. 
Viewing the charged black holes as near horizon limits of rotating brane configu- 
rations, the system is found to be metastable to one of the branes splitting from 
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the system of N coincident rotating branes. This is precisely the picture expected 
from weak coupling arguments where the locally stable vacuum can decay by one of 
the scalar eigenvalues tunnelling out into the unstable region. The qualitative phys- 
ical resemblance between weakly coupled gauge theory and dual gravity at strong 
coupling lends further support to the idea that some properties of gravity may be 
encoded in and extracted from weakly coupled gauge degrees of freedom. 



However, while pointing out the similarities between the weak and strong cou- 
pling phase diagrams, we should also emphasize the qualitative differences. These 
qualitative differences are the outcome of our analysis at low temperatures in the near 
critical region, and continue to pose intriguing questions. Reinterpreting the phase 
diagram of [26] in terms of grand canonical variables, /x and T, it has been pointed 



out in [IT, [34]] that the phase diagram at strong coupling has particularly interesting 



features at low temperatures and near the critical chemical potential. Specifically, 
only when all three chemical potentials are equal ^1 = ^2 = ^3 = the black hole 
instability line meets the first order transition line at T = and ^ = 1/R. At this 
point an extremal black hole solution with zero horizon radius becomes preferred over 
pure AdS. When the chemical potentials are unequal however, e.g. /i 2 = A*3 = 0, one 
finds that the black hole instability line and the first order line meet at in = 1/R 
with TR = 1/71*. At this point the charged black hole horizon shrinks to zero size 
and the actions of thermal AdS and charged black hole coincide. It is not known 
what happens to the gravity dual below this temperature. The picture we have found 
in weakly coupled field theory is quite different, since the instability line meets the 
li = 1/R line precisely at TR = in all cases. In all cases, at weak coupling, the 
instability line approaches the 11= 1/R line exponentially at low temperatures. 



7. Discussion and future questions 



In this paper we have investigated a corner of the phase diagram of M = 4 theory 
with R symmetry chemical potentials on S 3 . The region of interest corresponds to 
chemical potentials close to their critical values, beyond which the theory is known 
not to have a ground state. By explicitly computing an effective potential for the 
light degrees of freedom in this region, we showed that at the critical values for the 
chemical potentials and at zero temperature, flat directions open up. In particular, 
we find moduli spaces parametrized by N, 2N and 3iV bosonic degrees of freedom 
for one, two and three critical chemical potentials respectively. These correspond to 
| BPS, | th BPS and | th BPS sates of the M = 4 theory, respectively becoming light 
at critical values for different numbers of chemical potentials. The counting of zero 
energy modes at the origin of the moduli spaces is consistent with this interpretation. 
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Using the above picture of an effective potential, we calculate the same at small 
temperatures at critical values of the chemical potentials. The small positive thermal 
masses for the light degrees of freedom allow us to move away from critical chemical 
potentials, and indeed to exceed the critical values, resulting in metastable ground 
states. In all these situations, the theory is found to be in a deconfmed phase with 
the eigenvalues of the Polyakov loop collapsing to a point on the circle. These 
metastable plasma phases are the extension to low temperatures, of the metastable 



states at high temperature discovered in [111 j. Putting together the results at low 



and high temperatures, the resulting phase diagram of the theory at weak coupling is 
summarized in Figure 3. The remarkable resemblance to the strong coupling phase 
diagram in Figure 4 is clear. However, there are crucial differences as well. The 
main difference to emerge from our analysis is that the width of the metastable 
region goes to zero only at zero temperature (exponentially). This means that the 
instability line meets the first order deconfinement transition line at // = 1/R and 
T = for all patterns of chemical potentials. The strong coupling picture is quite 
different. The only situation where the first order Hawking-Page line meets the 
black hole instability line at T = 0, is when all three chemical potentials are equal. 
In all other cases, the instability line and the Hawking-Page line meet at a finite 
temperature. Despite these differences it is quite remarkable that the metastability 
discovered first in weakly coupled field theory, does appear to have a strong coupling 



analogue [24 



There are several related questions to pursue, along the lines of those addressed 
in this work. One obvious generalization of our results would be to the /3-deformed 
theory which has three global U(l) symmetries. The free theory will have the same 
phase structure as the M = 4 theory. However at non-zero weak coupling, the 
theories will differ. Due to lower supersymmetries, we expect that turning on different 
numbers of critical chemical potentials may result in scenarios quite different from 
the M = 4 case. Particulary interesting is the question of the existence of a stable 
Bose-Einstein condensed state at weak coupling. It may be possible to look for such 
a phase by turning on, say more than one critical chemical potential, followed by a 
small temperature. Looking for such an exotic phase in the dual gravity setup would 
be extremely interesting. 

In JT2|], it was pointed out that in the limit, T — > 0, /x — > R^ 1 with R J 1 _ fixed, 
the theory reduces to various quantum mechanical sectors, for e.g., the ferromagnetic 
Heisenberg spin chain. It is interesting to take this limit in our calculation of the 
grand canonical partition function and make contact with the results in recent studies 
of integrability in planar M = 4 SYM. 
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A. Scalar and Vector Laplacians on S 3 

Below we review some relevant aspects of the spherical harmonic decomposition of 
fields on S* 3 . To begin with, consider the kinetic terms for the gauge field: 

£ (gauge) _ rj r ip^pv = J d 4 xy/^Tr(-lA„{D 2 + A)A>*-±(D A + V i A i ) 2 ). 

(A.l) 

The Laplacian A = VjV* on S 3 depends on the tensorial nature of the object on 
which it acts. For example, on the vector component 

AA l = VjV j A l - RjAi , (A.2) 

where Rij is the Ricci tensor of S 3 . For the component A , which is a scalar on S 3 , 
the Ricci tensor part is not present and V 2 is equivalent to the scalar Laplacian. The 
eigenvectors of the scalar Laplacian are spherical harmonics Yn labelled by angular 
momentum quantum numbers, £ G Z > with 

AY, = -R- 2 £(£ + 2)Y e (A.3) 

and degeneracy (£ + l) 2 . The vector field Ai can be decomposed into the image 
and the kernel of : A 1 = B l + C % with ViB 1 = and C l = V 1 / '. These modes 
have different eigenvalues for the vector Laplacian. Firstly, those in the image of V 1 , 
i.e. C\ are given by V l Y t with i e Z > 0, which satisfy 

AV% = -R- 2 £(£ + 2)V% . (A.4) 

The remaining modes B\ in the kernel of V 1 , are spanned by V£, also labelled by 
the angular momentum £ G Z > 0, with 

AVj = -R- 2 {£ + l) 2 V; (A.5) 
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and degeneracy 2£(£+2). Finally, we will need the Laplacian on fermionic modes. For 
2-component real fermions 8 on S 3 , the fermionic Laplacian has eigenvalues R~ 2 (£+^), 
with £ e Z > and degeneracy £{£ +1). 

The gauge field modes Bi and Cj decouple from the chemical potential and at 
Gaussian order the relevant terms in the Lagrangian are 

^Tr B t [-bi- AW + v 2 ] B l + ^Tr B % [-D 2 - A«' + y? 2 ] B l , (A.6) 

where the superscript reminds us that the Laplacian is for divergenceless vectors and 
scalars, respectively. The prime means that the £ = mode is missing. When these 
fluctuations are integrated out their contribution to the effective potential is of the 
form 

j N oo 

-^ilogdett-^ + ^fe) 2 ] , (A.7) 
P ij=i i=t 

where di is the degeneracy of the modes with angular momentum quantum number 
i, so 2£(£ + 2) for B { and (£ + l) 2 for Q. The integer £ is the lower limit on the 
angular momentum quantum number, so 1 for both Bi and Cj. The quantity ee_ is 
the energy of the mode, so equal to 

B t : e e =^R~ 2 (£ + l) 2 + V 2 , 

A.8 

Q : e t = v /R~ 2 £(£ + 2) 2 + v 2 . 

The contributions from all the other modes have the same form for some set of 
{de,£o,ee}. 
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